Vibrational Enhancement of the Effective Donor - Acceptor Coupling 
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The paper deals with a simple three sites model for charge transfer phenomena in an one- 
dimensional donor (D) - bridge (B) - acceptor (A) system coupled with vibrational dynamics of 
the B site. It is found that in a certain range of parameters the vibrational coupling leads to an 
enhancement of the effective donor - acceptor electronic coupling as a result of the formation of 
the polaron on the B site. This enhancement of the charge transfer efficiency is maximum at the 
resonance, where the effective energy of the fluctuating B site coincides with the donor (acceptor) 
energy. 

PACS numbers: 34.70.+e, 73.40.Gk, 82.39.Jn 



Molecular electronics is progressing so rapidly that it 
is now possible to do measurements and assembly at the 
level of individual or few molecules [ij, |2| . Charge trans- 
port is known to occur in a wide range of linear chain 
molecules including the DNA double strand molecules. 
For DNA, it is believed that the charge transport phe- 
nomenon! is involved in the protection of the DNA en- 
coded information against the oxidative damage [3J- As 
the DNA molecule is essentially a dynamic structure on 
the time scale of charge transport, one expects that vi- 
brational dynamics to play an important role for DNA 
electronics, and, in general, for any property of biological 
molecules because biological functions of life are associ- 
ated with molecular motions but not to the static or dead 
structure (i.e., equilibrium positions of all atoms). 

In this paper we are interested in an one-dimensional 
DNA wire or bridge (B) connecting a donor (D) and an 
acceptor (A) sites. Usually, the bridge consists of N sites 
with one state per site (see the abundant literature de- 
voted to this issue in the Refs. p,IE I S[3,[5 l^ llolllll ll3- 
HQHEHEHIH ETSl E3) , and theoret- 
ical analysis of this problem requires to solve a system 
of N + 2 non- linear coupled equations. Unfortunately, 
such a problem cannot be solved analytically for N ^S> 1, 
and we have to recourse to numerical solution. However, 
many insights and essential features of the dynamics can 
already be gained and captured by studying a simple 
three sites: D - B - A. Generally speaking, the interac- 
tion between the donor and acceptor involves all states 
of the bridging subsystem. This bridge mediated inter- 
action can be characterized, under certain conditions, by 
a single energy dependent parameter - effective coupling 
- which plays the key role in the charge transfer. For 
small systems the phase coherence of charges is main- 
tained over the entire system, and the quantum effects 
are crucial in determining the system properties. On the 
contrast, for the long N ^> 1 bridges, fast relaxation 
processes result in a strong dephasing between charges 



in the system. Therefore, this leads to a rapid falloff of 
the off-diagonal elements of the density matrix such that 
the diagonal elements can described by a set of kinetic 
equations [ftI Eif . 

Our concern in this paper is to investigate a Hamil- 
tonian model describing the D - A coupling under in 
the presence of dynamic structural fluctuations. Such 
local fluctuations, including local vibrations, twist mo- 
tions, radial deformations and hydrogen-bond stretching 
or opening, are known to strongly influence charge trans- 
fer in DNA molecules [2^, |2(| |271 . For simplicity, we 
consider a three sites D - B - A system where electronic 
degree of freedom is coupled to an effective local vibra- 
tional degree of freedom. 

Let us assume that initially the charge is entirely lo- 
calized on the donor site with energy e. Then owing 
to the non-zero overlapping integrals of the electronic 
wave functions between the two neighboring sites, the 
tunneling of the charge takes place from the donor to 
the acceptor site with the same energy. Denoting by 
{|d>,|6>,|a>} the localized states on the the donor, 
bridge and acceptor, respectively, the Hamiltonian of the 
bridge - mediated charge transfer between the donor and 
acceptor acquires the form, 
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p mu> r , ,, 
H = H e +^ + — — + kr\b><b\, (1) 

in which the bare electronic part reads as, 

H c = e [\dxd\ + \a><a\]+e b \b><b\ (2) 
+v db [\dxb\ + \b><d\] 
+v ba [|6><o| + |o><6|] , 

where e is the one-site energy of the donor and acceptor, 
ef, the one-site energy of the bridge, m (m ~ 300amtt) is 
the mass of the bridge base pair, r is its radial displace- 
ment in the localized vibrational mode with frequency u>, 
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the momentum p is conjugated to r, and k is the electron- 
localized vibration mode coupling constant. The local- 
ized bridge mode can be treated classically since the cor- 
responding vibrational displacement amplitude is larger 
than zero point quan tum fluctuations for characteristic 
DNA parameters lUii (see also the following). 

The frequency of typical vibrations in DNA ranges lu ~ 
10 11 — 10 12 s _1 , and the scale of electronic overlap inte- 
grals between base pairs in DNA is v — y/2(v 2 b + v 2 a ) ~ 
0.2 eV leading to an electronic characteristic frequency, 
w//i~3x 10 14 s _1 . As a consequence of the small (adi- 
abatic) parameter, vu)/h <C 1, the slow vibrational and 
fast electronic motions can be decoupled. Therefore, to 
solve the problem of bridge-mediated charge transfer be- 
tween the donor and the acceptor, we employ the adi- 
abatic procedure to eliminate the slow vibrational mo- 
tions and derive an effective Hamiltonian for fast elec- 
tronic motions. To proceed, we take the wave function 
of the charge in the form, \^(t) >= a(t)\d> +c b (t)\b> 
+c a (t)\a > , where c n (t) are the time dependent ampli- 
tudes of the probability to have the charge at the nth 
site. From the Hamiltonian we arrive at the follow- 
ing equations of motion for the quantum amplitude c n (t), 




e v d b 
Vdb £b + kr v ba 

V ba £ 

and for the classical dynamic mode, 

rf2? " 2,1 IS 

m—rr = —moj r — k c& 
dt z 
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Next, we seek for stationary solutions of the form c n (i) = 
c n e~ lEt / h oscillating with frequency E/h. To work with 
dimensionless quantities, we use from now the dimension- 
less variables it, a and k defined in Table [I] Using c„(t) 
in the equations of motions, we find that the station- 
ary bridge displacement is r = — kc 2 and the stationary 
probability distribution is, 
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c 2 = 2u 2 /(l + 2u 2 ) 

^ = l/[(l + ^)(l + 2« 2 )] 



(5) 



where the eigen-energy u satisfies the characteristic equa- 
tion, 



4u 4 + 4(k 2 - cr)u 3 - 2<ru - 1 = . 



(6) 



Solutions of this equation provides the ground u g plus 
one or three (depending upon values of k and a) excited 
energies of the "polaron" , i.e. the state created by the 
charge coupling with the DNA structural deformation. 
These stationary polaron solutions allows us to eliminate 
from JJJ the structural deformation r and to obtain the 
effective Hamiltonian for the charge transfer as, 
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(fle-el) = A(u) \bxb\ 



TABLE I: Definition of dimensionless variables. 



definition 


vari 


able 


energy scale 


v — 




length scale 


c = 


^Jv/mui 2 


coupling asymmetry 


v = 


Vdb 1 Vba 


dimensionless polaron energy 


u — 


(E-s)/v 


dimensionless energy barrier 


a = 


(e b - s)/v 


reduced electron- vibration coupling 


k = 





+ 



i 



\dxb\ + \bxd\] 



[|&Xet| + |a><6|] , 



(7) 



where, A(w) = a — k 2 c 2 {u), is the renormalized effective 
energy of the bridge due to electron - vibration coupling, 
and 1 is the unit matrix. To calculate the effective charge 
coupling between the donor and the acceptor, we have 
to solve first the time dependent problem to determine 
the probability of charge transfer defined as, Pd^ a {t) = 
| < a|*(t) > | 2 , where > is the solution of the 

Schrodinger equation, 

dt 

with the initial condition, \^(0) >— \d >. It is easy to 
show that, when t-^oowe have Pd^ a {t) — kdat, where 
the charge transfer rate kda is given by the Fermi Golden 
rule, 



i 5(u 2 - u i) ; ° > k 2 cI 



kda = -^\ H da\ 2 X 



(8) 



d(u 3 -tii) ; cr < « 2 c 2 



with the (dimensionless) effective donor - acceptor cou- 
pling E3, 



\H d a\ 2 



4(1 + rff 
1/ [(ti 3 - tii)(ti 3 - u 2 )] ; cr > k 2 cI 

(9) 

l/[(u 1 -U 2 ){U3-U 2 )] ; cr < K 2 cl 

where Ui are the eigen-energies of H e f{ given by u\ = 
and 2it2.3 = A =p \/A 2 + 2. Finally, we end up with the 
effective D-A coupling given by, 
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2(1 + t? 2 ) 2 VA 2 +2 



[|A| + VA 2 + 2J ' 



(10) 
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and the ratio p, allowing to measure the effect of vibra- 
tions on the D-A coupling, reads: 



the distribution of charge density reduces to, 



p(u,a,n) 



\H d a(u,rj,a, k)\ 2 
\H da (u, V ,a,0W 



A| + VA 2 + 2 



-,1/2 



A 2 



(11) 



The expressions (|10|l and are the main results of 
this paper. They provide close formulas for evaluating 
how the dynamical disorder affects the effective donor- 
acceptor coupling in various situations. As a direct ap- 
plications of our main findings, we consider the following 
illustrative examples. 

Charge density versus polaron energy: 

It results from Eq.JSJ) that the charge densities on the 
donor, bridge and acceptor are even functions of the po- 
laron energy u. For all u, the ratio of the density of the 
donor to that of the acceptor is equal the square of the 
asymmetry energy 77 (see the Table HJ. For u = 0, there 
is no charge on the bridge site, and the charge density 
is distributed between the donor and acceptor sites in 
proportion of r\. In contrast, for the limits of very high 
(or low) polaron energy when |u| gets larger, the charge 
density decreases considerably on the donor and accep- 
tor sites while it gets higher on the bridge site leading 
hence to small charge transfer efficiency. These features 
are illustrated in Fig. ^ for the energy asymmetry pa- 
rameter -q — 0.5. Likewise, at the resonance u m where 
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FIG. 1: Stationary charge density in Eq.JSJ versus energy for 
the coupling asymmetry r] — 0.5. The quoted letters "d", "b", 
and "a" stand for the donor, bridge and acceptor, respectively. 

the renormalized effective energy of the bridge is equal 
to zero (the donor/acceptor energy), 



A(u m ) = 



u m = ± 



V^K 2 - &) ' 



(12) 



V 2 c 2 a (u m ) 
ct/k 2 



(13) 



c 2 (u m ) = [l- C 2 K0]/(l + r, 2 ) 



At this resonance point, the charge density on the bridge 
site decreases either upon approaching to the bare reso- 
nance for a — > or by increasing the electron-vibration 
coupling parameter k above \fa. 

Effective D-A coupling versus u: 

The effective D-A coupling scales as, \Hd a \ ~ 1/A 
for A 3> 1. However, at the resonance A = de- 
fined in Ea. lfT^jl . the \Hd a (u, r],a,n)\ 2 , and thus the ratio 
p(u, a, k), attain their maxima given by, 



\H da (u m ,a, k)Y 



4(1 + ?7 2 ) 2 



and, 



p(u m ,a,n) 



(a + 



(14) 



(15) 



Simple inspection of this equation shows two character- 
istic features. First, at the resonance the effective D-A 
coupling is enhanced by the coupling to structural dy- 
namics, and second, as it is illustrated in the Fig. |3 the 
enhancement factor p of the effective D-A coupling due to 
electron- vibration interactions increases with energy bar- 
rier a. To rationalize these observations in terms of the 
polaron energy and electron-vibration coupling, we have 
depicted in Fig. the ratio p(u, a, k) as a function of u 
for increasing values of a and n. Two different regimes 
can be distinguished: 

• k < ^fa: below the resonance value, the effective 
D-A coupling is a monotonic increasing function of 
the polaron energy \u\. 

• k > y/a: the effective D-A coupling increases for 
\u\ < \u m \, attains its maximum at the resonance 
M = |tt m |, and decreases for \u\ > \u m \. As a 
consequence of \Hd a (u, cr, 0)| ~ 1/cr, both the max- 
imum of p at the resonance and its limit at high |u| 
increase with a. 

Effective D-A coupling versus k: 

As we have discussed above and illustrated in Fig. [21 
the electron- vibration may lead to an increase or decrease 
of the effective D-A coupling depending on the value of 
|u| and the regime of a. Similarly, Fig. [3] displays the 
enhancement factor p at the polaron ground state as a 
function of the electron-vibration coupling k. It is clear 
that there is a certain threshold value K c (u g ) below which 
the electron-vibration coupling leads to enhancement of 
the effective D-A coupling, and above which the effective 
D-A coupling is drastically reduced affecting hence the 
charge transfer efficiency. 

As above, two different regimes can be distinguished: 
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• a < 1: Equation © has two distinct roots corre- 
sponding to the ground and excited states, respec- 
tively. The maximum enhancement, p(u rrll a, n c ) 
given by Eq. I|15|) , is attained at the resonance A = 
0, where u g = — V2/2 (and the excited state, — u g ) 
and K c {u g ) = \/2c7, obtained from the combination 
of Eqs.® and JEJ). In this regime, the ground 
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FIG. 2: Effective coupling ratio in Ea.ljllfl as a function of 
the energy. Filled circles correspond to p{u g ,a,K) at the 
ground state energy u g . The quoted numbers correspond to 
the electron-vibration coupling values, i.e., k = 0, 0.5, 1, 2, 5. 
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FIG. 3: Effective coupling ratio in Eq. lllH at the ground state 
energy u g as a function of the electron-vibration coupling. 
The quoted numbers correspond to the energy barrier values, 
i.e., a = 0.5, I, 2. 

state coincide with the resonant energy, u g = u m . 

• a > 1: There is an interval, K c (u g ) < k < n max , 
within which Eq. JBJ admits four distinct roots (the 
lowest one corresponding to the ground state) and 
out of which it has two distinct roots. In this case, 
the ground state is no longer resonant, u g ^ w m , 
but two excited states coincide with the resonant 
energies given by u m — ±^/2/2. As a result, the 
maximum enhancement p(u g ,a,K) < p(u rn , a, n c ) 
as A(u 9 ) 7^ 0. For instance, for a = 2, the 
four distinct roots interval is 1.9336 < k < 2.175 
with k c = 1.9336 and the maximum enhancement 
p(u g ,cr 7 K c ) — 3.198 is obtained for u g = —1.097. 
At the resonance A = for a = 2, we have k = 2, 
u g = —1.707, A(u g ) = —\[2 and p(u g ,a,n) — 
1.596. 

In summary, we have shown that the electronic cou- 
pling with the vibration dynamics of the bridge results 
in a formation of a polaron that may, under certain condi- 
tions, leads to an enhancement of the charge transfer effi- 
cency. Figures |21 and [3] show that the enhancement factor 
p is greater than one for a wide range of the energy barrier 
sigma and the electron-vibration coupling k. These find- 
ings are very suggestive for the issue of charge transport 
assisted by structural dynamical along the DNA chain. 
To study the basic mechanism of vibration enhancement 
of charge transport we have focused in this work on the 
simple three sites model with a single harmonic struc- 
tural dynamic mode (reaction coordinate). Meanwhile, 
the method employed in this work is not limited to this 
model and the extension of the theory to several sites and 
anharmonic reaction coordinates (see e.g., Refs-IUI^), 
and several resonance states can be handled within the 
framework developed in Ref. Nevertheless, further 

theoretical studies need to be conduct along the ideas 
outlined above in order to gain a better understanding 
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of charge transport properties in biological systems and 
technological applications of significant importance. 
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